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I. INTRODUCTION 


A. METAL MATRIX COMPOSITES(MMCs) 

Discontinuously reinforced Metal Matrix Composites(MMCs) represent a 
group of materials that combine the strength and hardness of the reinforcing phase 
with the ductility and toughness of the matrix. Aluminum alloys{Al) reinforced 
with Silicon Carbide(SiC) in particulate, platelet, or whisker form and fabricated by 


powder metallurgy methods are receiving a great deal of attention from researchers 


and engineers. 


B. THE ANALYSIS OF MMCs(Al/SICs) IN THE STABLE MEMBER 

Previous research has evaluated the dimensional stability, and the thermal 
and mechanical properties of several Al/SiC MMCs in stable—member applications 
for missile inertial guidance system. The results reveal that, although the candidate 
materials from a powder blend of SiC and Al alloy consolidated by Vacuum Hot 
Pressing{ VHP) into cylindrical billets, followed by Hot Isostatic Pressing( HIP) for 
full densification, have better isotropic properties, these MMCs show some micro 
creeping behavior in service. The microcreep of these MMCs will affect the 
dimensional stability of stable members, and is likely to arise from two sources: (1) 
creep of the metal matrix caused by internal stresses creep conditioned by externally 
applied stresses and/or (2) phase transformations during the creep condition. A first 
step toward understanding the cause of the dimensional stability problem is to 
analyze the influence of internal stresses. The internal stresses are thermal residual 


stresses when a Metal Matrix Composite is cooled down to room temperature from 


the fabrication or annealing temperature. These residual stresses are induced in the 
composite due to the mismatch of the thermal expansion coefficients between the 
aluminum alloy matrix and the silicon carbide particles. The model, based on 
Eshelby's for mismatch problems with simplified linear elastic material behavior for 
both particles and matrix, has been used to solve the problem of thermal residual 


stress [ref 1]. 


UC» ~PURP OSE 

The purpose of this thesis is to calculate the thermal residual stress of Al/SiC 
using theoretical methods and then to determine what influence this thermal 
residual stress has on creep deformation. In order to do this, the thermal residual 
stresses will be determined by focusing on elastoplastic matrix and elastic inclusions 
rather than on elastic matrix and elastic inclusions. Therefore, the case where both 
phases{ matrix and inclusions ) are perfectly elastic will be treated first. Then. the 
theoretical results will be compared with previously obtained experimental data for 


Al/SiC composites. 
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Il. LITERATURE REVIEW 

A. MMCs 

Metal matrix composites(MMCs), including eutectic composites, are becoming 
important in applications as structural components which are to be at intermediate 
and high temperatures. When MMCs are fabricated at high temperature or annealed 
at certain high temperatures, the MMCs have undesirable properties, such as low 
tensile yield and ultimate strengths. Those results are mainly due to residual 
stresses that are caused by the mismatch of the thermal expansion coefficient 
between the matrix and fiber. The residual stresses so induced have been observed 
in tungsten fiber/copper composites and in SiC whisker/6061 Al composites and 
were based on analyzed a one—dimensional model for a continuous fiber system or 
spherical particle system. The model. based on Eshelby's equivalent inclusion 


method, has been used to solve the problem of thermal residual stresses[ref 2]. 


B. ELASTIC MATRIX AND ELASTIC INCLUSIONS 

Let's consider a two—phase suspension of finite extent in which one phase 1s a 
matrix while the other is in the form of inclusions. It is assumed that the inclusions 
can be treated as identical ellipsoids with corresponding axes aligned. Further, the 
inclusions are assumed to be randomly distributed in the matrix in such a way that 
the suspension is homogeneous on a macroscopic scale. The task of the present 
section is to solve Eshelby's transformation problem for the case where the 
fractional volume f of the inclusions is finite. Following Eshelby [ref 2], we consider 
first an auxiliary problem. what is the resulting stress o. in an inclusion when every 


inclusion, having isotropic plastic moduli equal to those of the matrix, is subjected 


; Melt 
to a uniform transformation strain €~ ° 
1} 


impossible since no detailed information is given about the spatial distribution of 


To find the exact solution appears to be 


the inclusions. All we may find in this situation are approximate average fields. The 


average constrained strain ee in any inclusion may be written in the form 
1} 
ev aeC ye (1) 


Here, ¢ © is the constrained strain in a typical inclusion which would be produced if 
i] 


it alone is embedded in an infinite elastic body. This term is given. according to 


Eshelby[ref ], as 


C g T | 


Cu ee 
1j ykl 1 


LO 
— 


W here D ad is Eshelby's tensor which depends on the aspect ratios of the ellipsoidal 
1) 
inclusion and Poisson's ratio.(Here and in the following we use the usual summation 
il 
ae 
) 


convention: the range of subscripts is form 1 to 3) is a contribution due to all 


the surrounding inclusions and the presence of free boundary; this term is entirely 


elastic, and is assumed to be effectively constant in and around the typical 


C 


inclusion. Then, it will be a good approximation to regard Cc as the average 
elastic strain in the matrix for the contribution of any one inclusion to the average 
field; outside it is negligible when there are a great number of inclusion in the 
material. For a homogeneous elastic body in static equilibrium, a volume integral of 


any component of the elastic strain associated with an internal stress should vanish 


provided that the integration is taken throughout the body. Thus 


From (1), (2) and (3), we have 


aCe T a | 
. = (] ai ac (4) 


The average internal stress in any inclusion is given by 


je ( 
1 iki ok] k | 


on 
a, 


W here a is the elastic moduli of the matrix (and also of the inclusions in the 
present situation) [ref 3]. Where the inclusions have elastic moduli differing from 
those of the matrix, the same approach as used by Eshelby [ref 2] will be applicable. 
The essence of this approach is terms of equivalent inclusions. Let — be the 


1) 
elastic moduli of the inclusions, and €~ be auniform transformation strain of every 
] 


j 
inclusion. Then the transformation strain of the equivalent inclusions, eT can be 
1] 

determined by solving a set of equations 

~ eC ap ‘ me ke: 

C e - — = © — 

ar k ] cel ) 1) kl cy ee (6) 

with ( 4 ). The average stress in any inclusion is given by ( 5 ), of course. Consider 
a composite specimen which is free of stress at a references temperature Ty. Let the 


matrix and the inclusions have isotropic linear thermal expansion coefficients a and 


* 
a , respectively. The elastic state when the specimen is subjected to a uniform 


temperature change, ATo=T—To5, 1s described by the situation where the matrix 
and the inclusions undergo uniform transformation strains a@aT& and 
ij 
* 
a OT& respectively. Here 6 is the Kronecker delta. The resulting state of stress 
I 


J 1) 
can be analyzed by putting 


+ 
cl 


2 (a ees (7a 
ij ij 


++ 
The overall thermal expansion coefficients a can be expressed as a function of a, 


* 


a .f{ and the aspect ratios of the inclusions. 


MO ow 
ijkl? ~ ijkl 
C. ELASTOPLASTIC MATRIX AND ELASTIC INCLUSIONS 
When the temperature change is great, one should expect plastic deformation, 
either in the matrix and/or in the inclusions, which relaxes the internal stress 
accumulated. An important practical case, the case where the plastic deformation 
occurs only in the matrix, will be considered in the following. Further, for 
simplicity, we assume that the plastic strain is uniform throughout the matrix and 
the plastically nondeforming inclusions can be treated as identical spheroids. Then, 
the plastic strain in the matrix, en may be written in consideration of a symmetric 


mode of plastic deformation and the volume constant law, as: 
Pos ee 
on lO Hedy * 2%) (8) 


#* 
Where €, isa temperature—dependent parameter. In this case, the overall strain € 
1) 


of the composite is given, with use of Equation wz ician by 
1) 1) 


* ¥ 
é = aaTS +€ Pate t 
1j l)o1j 1] 


(2h) 


Where e! is calculated from (4) and (6) with 
ij 


a = a — ajaTs. —€ 2 (10) 
Now, the temperature dependence of €,, will be found through an energy—balance 
consideration. Since energy dissipative process is involved in plastic deformation, 
the decrease in the elastic energy, —éEe], corresponding to a virtual plastic 
deformation 6€,, should be balanced against the dissipation energy. If the matrix 


material has a constant flow stress o, in a uniaxial deformation, the condition is 


written as: 
—(EKe] = (1-f)o, | 6, | (0,,>0) (11) 
Pore Pye — Hor IT" 


* 
= Cres + Cae, (a — a)aT 
* 
+ C3(a _ a)*aT? ibe) 


It is noted that when the composite is elastoplatic the following condition should be 


satisfied: 


(Eel) > (Eel)ry 


—~J 


Where the subscripts I and II refer to the entirely elastic case and the elastoplastic 


case, respectively[ref 3}. 


Dee Feline 

If a stress is suddenly applied to pure metals, some solid solutions and most 
engineering alloys at a temperature near or greater than 0.5Tm(where Tm is the 
absolute melting point of the metal or alloy) deformation proceeds as shown in Fig. 
2.1. The initial application of stress causes an instantaneous elastic strain €9 to 
occur. If the stress is sufficiently high an initial plastic deformation €, also occurs. 
At low temperatures, significant deformation ceases after the initial application of 
stress and an increase in stress Is required to cause further deformation. At elevated 
temperatures, deformation under a constant applied load continues with time. The 
early stage of such deformation called primary creep is characterized by an initially 
high creep rate d€/dt which gradually decreases with time. Eventually a linear 
variation of creep strain accumulation with time is observed. This steady—state 
creep region is characterized by a constant minimum creep rate. Steady—state creep 
rates depend significantly on stress and temperature and are used frequently to 
compare the creep resistance among alloys. After significant deformation in 
steady—state creep, necking occurs or sufficient internal damage in the form of voids 
oI cavities accumulate to reduce the cross-sectional area resulting in an increase in 
stress and creep rate. The process accelerates rapidly and failure occurs. This region 
of the creep curve is called tertiary creep. Fig. 2.1 also shows the derivative of the 
creep curve or the creep rate curve. Fig. 2.1 represents the most common creep 
strain/time behavior and Fig. 2.2 illustrates the strain rate/time behavior at 


temperatures below approximately 0.3Tm. In Fig. 2, only transient creep is 


observed which is characterized by a continuously decreasing creep rate that 
approaches zero as the inverse of time. Such low-temperature creep behavior is 
called logarithmic creep. Fig. 2.2 also illustrates the creep rate/time behavior of 
alloys that exhibit a continuously increasing creep rate in the early stages of creep. 
Such alloys have been designated Class 1 alloys. Creep by viscous dislocation glide 
results when dislocations glide or move in slip planes under the action of an applied 
stress. These dislocations drag along solute atoms attracted to the strain fields of 
the dislocations. In order for the dislocations to move the solute atmospheres must 
diffuse in the direction of dislocation motion. At lower temperatures(~0.3Tm) or at 
higher stress levels, creep occurs by Syemelly, Sauk’ dislocation glide. Under 
these conditions dislocations can overcome barriers to motion without dislocation 
climb. The localized motion of short segments of small dislocations is important in 
overcoming barriers. Because of the larger contribution of stress in thermally 
activated dislocation glide a temperature and stress dependence for creep different 


from that for dislocation climb or viscous glide will be observed.[ref 6] 
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Fig. 2.1 Variation of creep and creep rate with time 
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Fig. 2.3 Typical creep curve 
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II]. THEORETICAL BACKGROUND AND DERIVATION 


A. THEORETICAL BACKGROUND 

To simulate Metal Matrix Composites(M{MCs), we consider an infinite 
isotropic body containing an ellipsoidal fiber as shown in Fig. 3.1. The infinite 
body(composite) experiences a temperature Tl — TQ. where T1 is assumed to be 
larger than the relaxation temperature of the matrix Tr and of course, Tr>TO. To 


make the problem simpler, the case of one fiber embedded in the infinite body is 


considered here. 





Fig. 3.1. Theoretical model 
Based on the model in Fig 2, let us consider the case of elastic 


inclusions—elastic matrix, elastoplastic matrix—elastic inclusions and creep 
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deformation: 


1. Elastic Inclusions and Matrix 
First let us consider the average internal constrain strain es To find 


E 


ond 


, we insert Equation(1) and Equation(3) into Equation(4) a find 


=(1-f)S wee} + ee. which represents the average internal constrain in any 
1) 1] 


inclusion in terms of transformation strain. Next, the average internal stress in any 


beet o 


inclusion can be represented by Equation(5). Here, es is an unknown term. 
Therefore, to find this unknown term we must consider Equation (6) and Equation 
(4). In Equation (6), * represents the inclusion, and ah represents the uniform 
transformation strain for every inclusion. Once we have determined 7 according 


to the Equation a fe we can determine the overall strain of the specimen, and 


am 
in accordance with Equation(7) we can determine the uniform transformation € 
i j 


o Elastoplastic Matrix and Elastic Inclusions 


The case for elastoplastic matrix—elastic inclusions is different from that 
of elastic inclusion and matrix. This is due to the fact that in this case plastic 
deformation occurs only inside the matrix, so we must consider the plastic strain in 


the matrix. Therefore we can find €P in accordance with Equation(S) and 
i 
* : 
Equation(10). From Ee = — Lele we can represent o =F(e and 
2 1) ij 
* 
et =F(¢ P). After representing them this way, 
ij 


1) 


Bel = — dig le™ 
Ms 1 j 
a2 
=C1€, 7+ C2€(a —ajaT + C2E(a —a) aT (14). 


— 
Cw 


Differentiating this with respect to € we can obtain the value of €, as shown 


below: 


RG yc ae 
One + ocr! y 


I 


Therefore, using the value of C to find the value € P, the value of o” can be 
ij i 
obtained, and inserting this value into Equation(16), we are ultimately able to 


obtain the total average stress in the matrix.[ref 9] 


facet | 
< a a Crk] ace Vp =—-f e ( 16} 
3. Creep 


Applying the value of the internal stress(above) to Equation(13). we can 
determine creep deformation. Following from this theoretical base, let us first 


calculate mathematically the thermal residual stress of Al/SiC MMCs. 


B. DERIVATION 
1. Elastic Inclusions and Elastic Matrix 
First let us calculate el for the case of elastic inclusions and elastic matrix. 
] 


In order to determine €  , we rewrite Equation(4) in matrix form: 
lj 


Ei 911.912 913.914 915 S16 |] Et Ey 
i - = 223 aa = a E> E> 
B37 qos 31 932 033 934 035 036 |] &3 €3 
ext =(1-4) 941.942 943 944 945 46 || E4 i 4 
€55 551 952 953 554 955 956 || Es E5 
| €66 N61 562 963 64065 O66 KL bo E66 | 
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Simplifying this matrix we obtain: 


en E 
€22 : 
Saar jeer ys 
ext =(1-f) ¢ 
Es € 
E66 E 


Further, simplifying this matnx we find: 


Ei) = €:f(1-£)Siy + €obSio( 1-f) + €3§S13(1-f) + fer} 

€55 = €1{(1-£)Sa1 + €2bSoo(1-f) + €3§S23(1-f) + feos 

€3§ = €1(1-f)Sgi + €ofSso(1-f) + €3$S39(1—f) + fesh 

C44 = €4§(1-f)S4q + fegh 

E55 = €5f(1-f)S55 + fess 

€66 = E6h(1-f)Se6 + f€6h (iq) 


Refer to Appendix (A) for the values of Sijx. 


Next, let us rewrite Equation(6) in matmx form. 


(I-v) vy vp 00 0 ei E; 
y (l-v) v 0 0 O €29 E2 
Gi vv (i-v) 0 0 O €33 |_| €3 = 
0 0 22 0 0 C44 E4 = 
0 0 O O24 O €55 ess 
0 0 0 O O 2p E66 Eg 
+ 
(l-v) v vy 0 0 O ei Eile 
y (l-v) v 0 0 0 E29 End» 
H]| vv (li-v) 0 0 O €33 |_| €33s 
0 0 0 24 0 0 C44 Eade 
0 0 0 0 Op 0 E55 €55% 
0 0 O O O 2B Ese} | € 


Simplifying the above matrix, we have the following: 


€71(1—v) J+-€35vJ+€33vJ=(1—-v) Gey} + Veob-+ vGesk—(1—v) Het 
_ ves — ec 

€j;vJ +€55(1—v) J+ €33vJ=vGeif + (1-v)Geb-+yGesh—rHe,} — 
(1-v) Hes — Hes}. 

€j;vJ+€39vI+€39( 1—v) J=vGe,}+ vGeés$+(1—v)Ges}—vHe,} — 
vHesS — (1-v)Hesh” 


d, 


E44 = 0 
é« = 0 
&e6 = 0 (18) 


From Equation (17) and (18), if we simplify and solve for ett we obtain: 
ij 


Kei}+Lesb+Mesh=(1-v) Het —v Hea —vHesh 
N€si+P€o5+ Qe;}=—vHe,} (1-v) Heo} —vHes} 
Reij+ S€o5+T€3h=—vHE ee sae 1—v)Hes} 


ef = Gh 
= est 
E 6k = ect ( 19 ) 


Now, we have 6 unknown values and 6 Equations. Therefore, applying Equation(19) 


to determine eT. 
i 


i i 
e f= —-YL — XM — HA (i+p) 
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W — VX 


€o$= 

ee ZU — XW 

ese ene 

€5§ = cu 

Eek = ect ( 20 ) 


Putting the above equations into simplified fomn: 


Es Y, 
E ae = E33 = X + 
ij E44 E 44% 
es} cant 
Es Es 


Here, refer to Appendix(C) for constants B through Z. Thus, inserting the value 
i 


€ into Equation(5), which is the average internal stress in any inclusion, we can 
ij 


determine the internal stress. 


On (i-v) vy vp 0 0 O]f BZ'+CY'+DX'-Z' 
0227 y (l-v) vy O 0 Of} C2Z'+BY'+DX'-Y' 
0331 |= G| vv (i-v) 0 O Of] FZ'+FY'+EX'-X' 
C447 0 0 O BH 0 0 0 
O55 0 OO 0 (0° 22350 0 
O66 0 0 0 0 O 2p 0 


2.  Elastoplastic Matrix and Elastic Inclusions 
Up to this point we have found the internal stress in the case of an 
elastic matrix and elastic inclusions. Now we will consider the case of elastoplastic 


matrix and elastic inclusions. Here also the method of determining the internal 


16 


and elastic inclusions. Here also the method of determining the internal stress is 
similar to the case of elastic matrix and inclusion. However in this case, since € P 
1) 


* 
the plastic strain, is included in the term el , the next case is more complicated. 
1] 


Ts 


Therefore, before we consider the relationship between € P and € * 
1) 1J 


, let us first look 
at = 
Bat - oo -€ Bl - 5 (5, 5 vr 6 5 I ( 6: i=j-1, i#j-0 ). Rewniting the 
1) 1 4) 1} 


Dove equation in iif ad we have 


E1fs b11 E11) 
C25 ae €22,, 
E344 3 €33 
=(a —a)aT — 
C444 (esa) Exar 
E55 % 5 €55 
et 566 Eger iy 


But according to Equation(8), e can be shown as follows: 


eur bi b14+401 doy 
C22), 512 619 +602 be 
€33, | =e 833 533+ 693 623 
E44, Pp 14 O1gt+04 doq 
C55), ( d15 Si5tdos dos 
Ege | A 616 Si6t+bo6 oe | 





From the above equation we can obtain the following: 


Ee 
-1/2€, 
eP. € 
ye p 
1j 0 
0 
0 (22) 


Therefore, inserting Equation(22) into Equation(21) and simplifying we obtain: 
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Ee * 
€ | =(a — a)aT + €)/2 


T* = 
= =(a — a)aT + €)/2 


+ + 
E e =(a — a)aT-—€p 


e "x0 (23) 


Next, using Equation(4) and Equation(6), if el is expressed as €]1, the result is the 
ij 
same as Equation(17). Inserting Equation(17) into Equation(6) and expressing el 
1j 
* 
as et we obtain: 
ij 


* % ' 

€ 7 =Qieit a ayank = Orese 
x x + 

Pie + Bren + eee 


* * * 
en + ee + meee 


E 


E 


=0 (24) 
Next, in order to find €jj, we insert the value of a into Equation(14) and simplify: 
ij 


en=Rie) tee Rs 
esRac\) epee) + Reeads 
econ Ieee Recon Beces 
€45=0 


€5=0 
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€56=0 (25) 


Therefore ort can be found in accordance with Equation(5). (The value of this ra is 
different from the value of oi for elastic inclusion and elastic matrix, the reason 
being that oi} for elastic inclusion and elastic matrix does not contain e.). Here, in 
order to determine the value of oi! we must first find the value of €.: From 
Equation(10), if we differentiate E,; with respect to € we can obtain Equation(11). 
Therefore, from 6Ee) =| 2Ci€,, + Co( a = a) aT] de, and —éEe] = (1-f)o, |e, 1; we 
can obtain the value of e Again, inserting this value of €, into Equation(22) we 
can find the value of eo Finally, by inserting the value of ous into equation(27) we 


obtain the value of iy. Utilizing Equation(10) again to find the value of €, 


* 
ving |f ens Re UeeD cera Alsace) 2 
0227 €95% X3A! + X4€p A' + €,/2 
= 2 744) €4 x 2 0 0 
055] E55 + 0 0 
766 |} E66 | | 0 0 


=— , [e5( ae — Xe) + €)A( mH Xo+e Ut Xi-XstXe ) 
+ A'*( Xi+X3+Xsz )] ( 26 ) 


Therefore, differentiating E,] in respect to €,: 
bE e} =—f( C2Al + 2EpCy )5, ( 2¢ ) 
And inserting Equation(27) into Equation(11) and simplifying obtains: 


f( CoA’ + 2€pC1 dE =(1-) oy ( 28 ) 
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If we find the value of €, from Equation(28). 


Seo ee 
€p=— 90,(¢ — aal + a {1 hoy (27) 
( 4: 6€,<0 ~ Cooling, — be, >0 ~ Heating ) 


If this value of €, is inserted into Equation(24) and (25), the value of an and the 
1) 
value of €~ can be found, and if these values are inserted into Equation(5) the value 
ij 
of oe can be determined. Therefore, inserting this value into Equation(16) we are 
1 


ultimately able to obtain total average stress in the matrix. Refer to Appendix( C ) 


for the constants and actual values. 


3. Creep 


Up to this point we have been determining the internal stress both in the 
case of elastic matrix and elastic inclusion and in the case of elastoplastic matrix 
and elastic inclusion. But in these two cases the actual effect on creep deformation is 
the internal stress in the case of elastoplastic matrix and elastic inclusion. 
Accordingly, by inserting the value of the internal stress into Equation (13), we can 
plot the microcreep deformation phenomenon as a function of time as shown in Fig. 


3.2. Refer to Chapter IV part A for more detailed information. 
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FIG. 3.2 MicroCreep as a function of time t 
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IV. RESULTS AND CONCLUSIONS 
Ao ee 
The thermo—mechanical data of the matrix and fiber for the theoretical 
calculations are obrtained from the [ref 13}. 


Annealed 2024 Al matrix: 


Emi eka 

a 47.5 MPa 

y =0.3 

a Se l0mo ian (31) 
SiC fiber: 

Biate7 Gea 

i= 02 


+ 
oe Lego) kK 
/d= 1,5 (32) 
Where the average value of the fiber aspect ratio(l/d) was used [ref 13]. The 


temperature drop AT is define as 
AT = T1—To (33) 


Where Tl is taken as the temperature below which dislocation generation 1s 
minimal during the cooling process and To is the room temperature. Thus, for the 
present composite system AT is set equal to —200K. From the data given by 
Equation (31), (32) and the use of < o a so, we have computed the stresses. 


Next, the thermal residual stresses, averaged in the matrix of SiC fiber/2024 Al, are 


to 
to 


predicted by (17) and the result on < o >, as plotted in Fig.6, where 33 denote 
the component along the longitudinal direction(z). The average theoretical thermal 
residual stress is predicted to be tensile in nature, and the average residual stress in 
the longitudinal direction to be larger than the average residual stress in the 
transverse direction [Table 5]. The fiber aspect ratio(a=1/d) of SiC fibers has been 
observed to be variable [Appendix A,B]. In the present mode] we have used the 


value of I/d, 1.5 to predict the thermal residual stress of the composite. 


TABLE 1: The value of € P (a=i,1.5) 
1) 


a(1/d) € P eP e P 
1} oo ao 

1 0.0017 0.0017 0.0017 

1.5 0.0014 0.0014 0.0028 


TABLE 2: The value of _ (a=1, 1.5) 





TABLE 3: The value of ej (cal Ss) 









<o 


<o >(MPa) 


55.245 095.245 99.245 
30.905 70.998 





<o, >| MPa) > MPa) 






TABLE 7: Comparison value of < a Sp with TODAS 
1 


a{|/d) OURe MODELS 'S ewe 
< 033 >m( MPa) < ©3233 >n(MPa) 
1.8 7 eee 5 6 67.894 


As seen in tables 1. 2, 3, 4, 5, 6, and 7 we can observe that as a increases, the values 


ee), © iy, €ij\, a4 La and < a > also increase. 


Also in Fig. 4.1 we observe that when the aspect ratio(a) increases, the value of the 


internal stress also increases with it. 


22 MPa 


70 
68 
66 
64 


62 
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' Fig. 4.1 Internal stress vs aspect ratio 
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And in Fig. 4.2 we can also see that the value of the interna] stress increases with an 


increase in temperature. 


co co cE I CH 





JOO 400 


foo 150 200 290 300 
ae Temperature change(deita 7) 


Fig. 4.2 Internal stress vs temperature change(T ) 


— 


_—— 


2¢ 


In Fig. 4.3, in the case of an ellipsoid (a=1.5), the values of internal stresses 


increase with an increase in volume fraction, 


eooa T= (th 


rT 





8. a 0,12 0,14 0. 16 0.18 0, 


Volume fraction 


Fig. 4.3 Internal stress vs volume fraction(f) 
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Next, if we compare the theoretically predicted value of the internal stress with the 
value determined by R. J. ARSENAULT and M. TAYA, we see that our model's 
value is larger than the value which is obtained by R. J. ARSENAULT and M. 
TAYA, as shown in Table 7. They used the material properties as follows: 
Annealed 6061 Al matrix: 

E=47.5 GPa 

ao—aio MPa 

y =0.33 

a =23.6x10°/K 
SiC whisker: 

Ej;=427GPa 

Ok 

ems i 

1/d=1.8 
Finally, Fig. 3.2 is a graph showing the strain as a function of time, when the aspect 
ratio is 1.5 and the value of the thermal residual stress is inserted into 
Equation(13). If we analyze Fig. 2.3 we can see that the creep deformation 


phenomenon is due to the internal thermal residual stress of the MMCs. 


B. CONCLUSIONS AND RECOMMENDATIONS 

The object of this research is to obtain the value of the thermal residual stress 
of an Al/SiC composite using Eshelby's theoretical model, and then to determine 
what effect this thermal residual stress value produces on creep deformation. 
Because exact creeping behavior of this material is difficult to determine when 


analyzing creep deformation, Andrade's model was used with properly selected 
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constant values. But through this theoretical approach the following conclusions can 


be obtained. 


1. Thermal residual stress due to the difference in thermal expansion 
coefficients may be estimated. 


a. By means of the theoretical model, we can see that the thermal 
residual stresses increase when the volume fraction of the inclusions 
increases. 


b. In the case of sphere inclusions, that is, when the aspect ratio a=1, the 
lateral stress is equal to the longitudinal stress. 


c. In the case of an ellipsoid inclusions where a=1.5 the longitudinal stress 
is greater than the lateral stress. 


d. The thermal residual stresses of Al/SiC composites increased when 
the value of the aspect ratio of the inclusions increases. 


2. Microcreep deformation can be estimated in the model by using the 
thermal residual stresses and the Andrade's model of creep 
deformation 


a. Microcreep deformation is due to the internal thermal residual 
stress of the MMCs. 


b. Dimensional stability of components will be influenced by the 
behavior of the composite. 


3 Recommendations 


a. Presently, this thesis has only dealt with average thermal residual 
stress from a overall point of view. and we need more detailed 
local residual stresses surrounding SiC particles should be 
estimated for analyzing creeping Sects: 


b. In the future, research should be done concerning relaxation during 
cooling. 


c. The creep behavior used in the current model is the Andrade's 
approximation. The real creep deformation for this Al/SiC 
composite should be further studied in order to get better 
understanding of the microcreep of the material based on 
approximation values. 
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APPENDIX A 
PROGRAM FOR VALUE OF S,,(a=1.5 — 5) 


1. PROGRAM(FORTAN) 


300 


90 


95 


REAL*8 NU,ALP,G, ‘mp, §1111,$3333,81122, $1133, $3311 
NU=0.3 2 ; 

WRITE (6,90). zz 
WRITE(6,95) ; 
ALP#1.5 
TMP=ZALP**2-1. ) 
G= (ALP* (ALP*DSQRT (TMP) -DLOG (ALP+DSQRT (TMP) )))/ (THP** (3. Ja.)) 
S$1111=#(3.*ALP**2) /(8.*(1.-NU) *TMP) 3 3 
++(1.-2.*NU-9./(4.*TMP) ) *G/(4.* (1-NU) ) | 
$3333=((1.-2.*NU+(3. *ALP**2-1. )/TUP-(1.-2. *NU+3. *ALP**2/TMP) 
+*G))/(2.%(1.-NU) ) 

S$1122=(ALP**2/ (2. *TMP) —(1.—-2*NU) -3*G/ (4. aTHP))/ (4. #(1.-NU) ) 
S$1133=-ALP**2/(2.*(1.-NU) *TMP) + (3. *ALP**2/TMP-(1.-2.*NU) ) *G 
+/(4.%(1.—NU) ) ! 

$3311=-(1.-2. *NU+1. /(TMP))/(2. a(1, -NU))+(2.~2. *NU+3./(2.*TMP) ) 
+*G/(2.%(1.-NU)) | 

WRITE (6,100) ALP, $1111 ,$83333,81122, S1133, $3311 


ALP#ALP+.1 | 
IF(ALP.GT.5.) GOTO 400 

GOTO 300 Ae 

FORMAT (1X, ' ALP $1111 $3333 ' $1122 $1133 
+ $3311') . Z 


FORMAT ( 1X, e SSS SS S655 SS SS SS SS SS SSS SSS SSS SSS SS SSS SS SSS SSS 


+uuuus ° ) 


100 FORMAT (1X,6(F10. 5)) 
400 STOP , 


END — , / 
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2S 
Dy Sik 


ALP 


$1111 


1 VALUE WITH ASPECT RATIO 


$3333 


$1122 


$1133. 


$3311 


1.50000 
1.60000 
1.70000 
1.80000 
1.90000 
2.00000 
2.10000 
2.20000 
2.30000 
2.40000 
2.50000 
2.60000 
2.70000 
2.80000 
2.90000 
3.00000 
3.10000 
3.20000 
3.30000 


3.40000 - 


3.50000 
3.60000 
3.70000 
3.80000 
3.90000 
4.00000 
4.10000 
4.20000 
4.30000 
4.40000 
4.50000 
4.60000 
* 4.70000 
4.80000 
*4.90000 


- 58078 
-58828 
~59495 
~-60089 
-60620 
-61098 
-61529 
-61920 
©62274 
-62598 
-62893 
-63164 
-63412 
63641 
63852 
~-64048 
-64229 
~ 64397 
64554 
- 64700 
-64836 
- 64963 
-65082 
-65194 
-65299 
-65398 
~-65491 


~65579° 


- 65662 
© 65740 
~-65814 
~-65884 
-65951 


-66014. 


~-66074 


39378 
37371 
-35518 
33804 
- 32215 
30739 
. 29366 
- 28087 
.26893 
° 25777 
24732 
- 23752 
- 22832 
- 21967 
221152 
20384 
19659 
18974 
18326 
17713 


eager 


~16579 
~16054 
°15555 
~15081 
~ 14629 
-14197 
- 13786 
© 13393 
~13018 
~ 12659 
°12316 
- 11987 
-11671 
-11369 


.01421 


~01603 
- -01764 


~01906 
- .02032 
°02145 
~02246 
©02337 
~02419 
002493 
-02560 


.02622 — 


-02678 
002729 
-02776 
-02819 
.-02859 
~02896 
-902930 
~02962 
02991 
~03019 
03044 
-03068 
-03090 
-O03111 
©03130 
-903149 
~03166 
-03182 
~ 03197 
°03212 
~03225 
003238 
~03250 


e 


08396 
08992 
.09549 
10071 
. 10560 
.11018 
.11448 
2118521 
212231 
.12588 
12924 
213241 
.13540 
.13822 
-14090 
014342 
.14582 
.14809 
.15025 
.15230 
.15425 
-15610 


315787. 


.15955 
.16115 
.16269 
.16415 
.16555 


-» 16689 


- 16818 
~ 16940 
017058 
017171 
~17280 


~01945 
-01623 
~ 01350 
©01117 
-00918 
~00747 
- 00600 
~00472 
-00362 
~-00266 
-00183 
~00111 
-00048 
-.00007 
-.00056 
-.00098 
-.00134 
-.00167 
-.00195 
-.00219 
-.00241 
-.00259 
-.00276 
-.00290 
-.00302 
-.00312 
-.00321 
-.00329 
-.00335 
-.00340 
-.00345 
-.00348 
-.00351 
-.00353 
-.00355 


APPENDIX B 


1. PROGRAM(MATLAB) 


*% 

SSPHERE (ALPHA=1.0) 

%& 
for i=#1:6; 
nu#i/10.; 
den=15.*(1-nu); ae 
s11(i)=#=(7-5*nu) /den; 
s12(i)=(5*nu-1) /den; 
813(i)=s12(i); : 
831(i)=s12(i):; 
833(i)=s11(i); 


é 
aos ——s ow Oem le lle 


& | 5 : 
XALPHA=INPINITY v 
& : 
den2=8*(i-nu); . ‘ 


‘$s11(i)=(5-4.*nu) /den2; 
8812 (i)=-(1-4*nu) /den2; 
ss13(i)=nu/ (2. *(1-nu) ); 
8831(1)=0; 
$s33(1)=0; 
end; 
resulti=(s11;s12;s813;s31; 833; ye 
result2=(ssll1;38s12;:ss13;ss31;8ss33] ' 
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2 Sik VALUE 









C2 €O Oo @ @® 
ee ef + o& oo 


Cc? 


SS © © © GC ©o 
- ° ° ° ° ° 
OH View wre 


(VY Gt ee W IN 


resultl = 


0.4815 
0.5000 
0.5238 
0.5556 
0.6000 
0.6667 


result2 = 


0.6389 
0.6563 
0.6786 
0.7083 
0.7500 
0.8126 


-0.0370 


0 
0.0476 
O22 er 
0.2000 
0-333 


-0.0833 
“0.0312 


0.0357 
0.1250 
0.2500 
0.4375 


-0.0370 


34 


0 
0.0476 
Oo. 12a 
0.2000 
O-3335 


0.0556 
0.1250 
0.2143 
0.3323 
0.5000 
0.7500 


-0.0370 


0 
0.0476 
0.1111 
0.2000 
0.3333 


o0O00 0 0 


0.4815 
0.5000 
0.5238 
0.5556 
0.6000 
0.6667 


oOoO00 0 OC 


| 
! 
, 
. 
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APPENDIX C 


PROGRAM FOR VALUE OF THERMAL RESIDUAL STRESS, 


P TF CC ol ~| 
Cp Gi Fie” ae a 


i. PROGRAM(MATLAB) 


xed 
&This program computes stress values and constants 


x 

% Initial values re 
x | ; 
nus0.3; 

em=47.5e93 

ei#é27 .0e9? 

sy=47 .5e6; 

£=0.2; 

dt=-200 

ap=23.6e-6; | 

as#4.3e-6; 

muz0.; 

% : 
% Eshelby's tensor components 

% 

g11=0.6786; 

s83320.0000; 

s12=#0.0357; 

$13#0.2143; 

s3120.0000; 

% . 

& Constants ¢ 
x . 
app=(as-ap)*dt; 

gzen/((1. +nu) *(1.-2. *nu))3 - 
h=eei/((1.¢nu) *(1.-2.*nu) ); 
j=(em-ei)/((1.+nu)* (1.-2.enu)) 5 
a=f+s11%*(1.-f£); 

bes12*(1.-£);3 

cusli3*(1.-f); | 
e=f+s33*(1.-f)3 

d=s31*(1.-f); 


t 


continued 


tleri+r2+tr3-ql+q2+q3; . 

t2=(ri¢tr2- SRE O/T a 

t3urd+r5+r6-pi-p5-p6; _ 

td=(ré+r5-2*r6) /2.-(pd+p5-2*p6) /2. ; ' 
t5=ar7+r8+r9-pl-p2-p3; . 

t6=z(r7+r8-2*r9) /2.-(plt+p2- 2*p3)/2.; 

X1l=g* ((1L.—-nu) *tl+nu® (t3+t5));: 

x2ug*((1.—nu) *t2+nu* (té+té6)): “- 

x3eg* ((1l.—-nu) *t3+nu* (t1+t5) ); 

x4mug*((1.-nu) *t4é+nu*® (t2+t6)): 

x5ag*((1.—-nu) *t5tnu* (t1+t3)):;: 

x6eg*((1.—nu) *té6+tnu* (t2+t4) ); 

Cl=((x2+x4) /2.-x6)*(£/2.); 

C2=((x1+x3) /2.4+x2+x4-x5+x6)*(£/2.):3 

% Temprature dependent. parameter 

% ae 

ep=-c2*app/ (2.%c1)+(1.-£) *sy/(2.*%¢1); 

% } , . 

% Elastic moduli of the matricies 

% : | . 

aaa=g*({(1l-nu) nu nu 0 0:0;nu (1-nu) nu 0 O O;nu nu (1-nu) 0 0 0; 

000 2.*mu 0 0;0 0 00 2.*mu 0:0 0 000 2.*mu)'; 

% ; >. 

% Plastic strain in the matrix 

% 

eljp=(-ep/2 -ep/2 ep 0 0 0)’ 

% | 

% Uniform transformation strain of inclusion 


% 
eijts=[apptep/2 apptep/2 app-ep 0 0 0)" 


é a 
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continued 


% ; 
eijt=(ql* (apptep/2)-q2* (apptep/2)-q3*(app-ep) .. 
p4* (apptep/2)+p5* (apptep/2)+p6*(app-ep) .. 
pl* (apptep/2) +tp2* (apptep/2)+p3* (app-ep) 0 0 0) 
eij_c=#(rl1* (apptep/2) +r2* (apptep/2)+r3*(app-ep) .. 
r4* (apptep/2)+r5* (apptep/2)+r6*(app-ep) .. 
r7*(apptep/2)+tr8* (apptep/2)+r9*(app-ep) O 0 Q) 


oo 


Average internal stress in inclusion 
ij_isaaa* (eij_c-eijt) 


Average stress in the matrix 


eereeu ee 


sijm=-f*sij_i 
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B=S,,(1-f) + f 

aS et 

F=S9)(1—-f) 
H=Ej/(1+v)(1—2v) 

K= aR ele ae 
M=J((1-v)D+vD4+vE)—vG 
P=J((1-v)B+vC+vF){(1-1)G 
R=J((1-v)F+vB+vC)—vG 
Y=J((1-v)B+2vD)41-1)G 
V=(QK —NM)/K 

X=(SK —RL)/K) 

Z=A H(1+v)(R-K)/K 

Y =(W-X V)/U 


2. CONSTANTS OF ELASTIC MATRIX AND ELASTIC INCLUSIONS 


C=Syo(1-f) + f 
E=S39(1-f) + f 
G=Ba/(1+v)1-27) 
J=(Em—E;)/(1+v)(1-2v) 
L=J(i—v)C+ vB 4rE ane 
N=J((1-v)D+vD4+vE)—vG 
Q=J((1-v)D+vD+vE)—vG 


S=J((l-v)F+vC+vB)—vG 
U=(PK Sik 


W=A H(141)\(N-K)/K 
V=Giiee MRK 

X =(ZU-XW)/(UY-XV) 

7, (-\ JX NEHA Clee 
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APPENDIX D 
FORMULA OF Sijx) FOR A VARIETY OF INCLUSION 








‘ 5 eye cre x ] f - 9 | 
Sin = S222 * : mm 1 2 — oo) B 
Pema lil 1) 4(1 — va S Genin 


5 a (Satz — 1) { Be Sort | | 
Ak: BRCUOMS Temas (bo oe & 
Be i. yn) ara ere ee) 








rea 3 
ates me Mee) 
; oa _ ! te ] Barr - 
o>R23 *) 2233 (ex = 1) ee 


2(1 — v9) (arr — 21) arr — WW) 


- (= 2a 





a 
. “5 = ‘ SS «= oof = ) 
Sians = Sag22 ii Sen t 2vq + oe -| 
2(1 — 0%) 21> = 1] 


where vy is Poisson's ratio of a matrix, cr is aspect ratio of a fiber 
(=//d), and g 1s given by 


a . 
Bo [arttar*® — 1)! -— cosh 
(at = 1) 2 € 
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—_ ¢ — oy 

1111 a gkee 

_ wl 

1122 151 = v) 

 &-!l 

1133 ear | 

_ w—l 

9211 [5(1I — v) 

— tf -— OV 

> ee — Vv 
5yvy — l 

9233 «44 a 

5y — 1 

~33iy dc — v 
5y -— 1 


40) 


_ o— 4v 
itt nd 
4y — 1 
1122 — Vv 

= Vv 
133. 21 — v) 
C aa 4vy — |] 
oot) Ol — VY) 
_ o5—-— 4 

S = 

9009 Ol — V 
9993. 4ll — V 


3322 


3333 
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APPENDIX & 
PROGRAM FOR CREE! DEFORMATION 


PROGRAM(MATLAB) 


epd=0.001; 

beta=0.0465; 

k=7.2e-7; 

for tte=ic9l; 

KV( CE) =e 

c=1,0*820>( (1Ct-1)7 307) 
eep(tt)=ep0*(1.tbetra*t~ (1/3) ) *exp(k*t) ;. 
end; 

plot(kv,eep) 

semilogx:; 

xlabel(’'Time t'),ylabel('*Strain eps’) 


10. 
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